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Abstract. High angular resolution Australia Telescope Compact Array (ATCA) observations of SiO 'thermal' 
millimetre line emission towards the two oxygen-rich, low mass loss rate AGB stars R Dor and Pup are 
presented. In both cases the emission is resolved with an overall spherical symmetry. Detailed radiative transfer 
' modelling of the SiO line emission has been performed, and the comparison between observations and models are 

, conducted in the visibility plane, maximizing the sensitivity. The excitation analysis suggests that the abundance 

' of SiO is as high as 4 x 10~^ in the inner part of the wind, close to the predicted values from stellar atmosphere 

' models. Beyond a radius of « 1 x 10^^ cm the SiO abundance is significantly lower, about 3 x 10~®, until it 

. decreases strongly at a radius of about 3 x 10^^ cm. This is consistent with a scenario where SiO first freezes out 

' onto dust grains, and then eventually becomes photodissociated by the interstellar UV-radiation field. In these 

(-H ^ low expansion velocity sources the turbulent broadening of the lines plays an important role in the line formation. 

^ Micro-turbulent velocity widths in the range 1.1 — 1.5 kms~^ result in a very good reproduction of the observed 

line shapes even if the gas expansion velocity is kept constant. This, combined with the fact that the SiO and 
CO lines are well fitted using the same gas expansion velocity (to within 5 — 10%), suggest that the envelope 
acceleration occurs close to the stellar photosphere, within < 20 — 30 stellar radii. 
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1- Introduction tion on this region comes from observations of infrared ro- 

^ , vibrational molecular lines in ab sorption, main ly tow ards 

The intense winds that low to intermediate mass stars de- the carbon star IRC+ 10216 re.g.. lKeadv RidgwavJ flB 

velop m their final evolutionary stage, as they ascend the [winters et ai. I Effln^ . infra red continuum emission from 

asymptotic giant branch (AGB), return enriched stellar circumstellar dust (e.g.. iDanchi et"Zlll99l . and SiO 

material to the interstellar medium. In addition to signif- niase r radio line em ission from excited vibrational states 

icantly contributing to the chemical evolution of galaxies, g [Cotton et al lEoO^ 
the mass loss will dictate the time scale for the future 

evolution of the star towards the planetary nebula phase. Information can also be gained from radio observa- 
Considering its importance little is known with certainty tions of suitable circumstellar molecular species towards 
about the mechanism(s) behind the mass loss. Partly this stars having low to intermediate mass loss rates. A major 
is due to a lack of observational constraints, in partic- survey of CO radio line emission from circumstellar en- 
ular close to the stellar photosphere where the wind is velopes (CSEs) around oxygen-ri ch AGB star s of different 
accelerated, and partly due to the complexity of the phys- variab ility types were done by [Kerschbaum fcO Jofsso^ 
ical/chemical processes involved. Among other things, the U2M) andjOlofsson et alj l|200^. These data were mod- 
formation of dust grains is thought to play an important elled in detail to derive stellar mass loss rates an d 
role in these radiatively driven winds. Most of the informa- terminal gas expansion velocities (|01ofsson et^|200^. 

Subsequently, a survey of 'thermal' SiO radio line emis- 

Send offprint requests to: F. L. Schoier sion, meaning emission from the ground vibrational state 

e-mail: fredrik@astro.su.se rotational lines which are normally not (strongly) maser- 
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ing, was done and the data were interpreted using a de- 
tailed numerical radiative traii sfer modelling presented in 
iGonzalez Dekado et all l)2003tl . 

An immediate conclusion from these large surveys is 
that the SiO and CO radio line profiles are different from 
each other. The SiO line profiles are narrower in the sense 
that the main fraction of the emission comes from a ve- 
locity range smaller by about 10 — 20% than twice the 
expansion velocity determined from the CO data. On the 
other hand, the SiO line profiles have weak wings, such 
that the total velocity width of its emission is very similar 
to that of the CO emission. Furthermore, it appears that 
the SiO line profiles change character with the mass loss 
rate, at low mass loss rates they are narrow with weak 
extended wings, while at high mass loss rates they be- 
come distinctly triangular. These 'peculiar' SiO line pro- 
files have been interpreted as being due to the influence 
of gas acceleration in t he region which produ ces most of 
the SiO line em ission llBuiarr abal et alTllQS B^. However, 
as illustrated in IConzalez Delgado et all l|2003'). the SiO 
lines are usually strongly self-absorbed also for low mass 
loss rate objects and this produces narrower lines. There is 
some controversy in the literature over the scale length of 
the acceleration region. In the only previously published 
works where thermal SiO emission have been observed to- 
wards AGB stars using interferometry, ISahai fc Biegind 
(|1993') find n o need for th e slow ly varying velocity fields 
introduced bv lLucas et al 

Other conclusions come from the detailed radiative 
transfer modelling of the SiO line data. This is in many 
respects a more difficult enterprise than the CO line mod- 
elling. The SiO line emission predominantly comes from a 
region closer to the star than does the CO line emission, 
and this is a region where we have fewer observational 
constraints. The SiO excitation is also normally far from 
thermal equilibrium with the gas kinetic temperature, and 
radiative excitation plays a major role. Finally, there ex- 
ists no detailed chemical model for calculating the ra- 
dial S iO abundance distribution. iGonzalez Delgado et alJ 
l|2003(l adopted the assumption that the gas-phase SiO 
abundance stays high only very close to the star, since fur- 
ther out the SiO molecules are adsorbed onto the grains. 
Beyond this the abundance stays low until the molecules 
are eventually dissociated by the interstellar UV radia- 
tion. This photodissociation radius, which is crucial to 
the modelling, was estimated using both Si O multi-line 
modelling and existing inter ferometer data ((Lucas et alJ 
ll992H Sa,ha,i Bieging but only for a few sources. 

The result of the radiative transfer modelling is a circum- 
stellar SiO abundance that is roughly the same as that 
obtained from stellar atmosphere equilibrium chemistry 
for low mass loss rate objects, and which declines with 
mass loss rate reaching an abundance about two orders of 
magnitude lower for high mass loss rate objects. 

Thus, there arc strong indications that 'thermal' SiO 
radio line emission is a useful probe of the formation and 
evolution of dust grains in a CSE, as well as of its dy- 
namics. Hence, circumstellar SiO line emission potentially 



carries information on the properties of the region where 
the mass loss of AGB stars is initiated. 

Presented here are results from high spatial resolu- 
tion imaging of the two M-type semiregular AGB stars 
R Dor and Pup. Both sources studied have low mass 
loss rates (M < 1 x 10^^ Mq yr~^) and low terminal ve- 
locities of their winds (< 6 kms~^). It has been suggested 
from hydrodynamical calculations (JVintcrs ct al. 2 0oT)bl 
120021 boO.^h that the main driving mechanism behind these 
tenuous winds is stellar pulsation and that dust plays only 
a secondary role. 

The observations, performed with the Australia 
Telescope Compact Array ^ (ATCA), are presented in 
Sect. 2. ATCA is an array of six 22 m dishes operating 
from 1.4 — 26 GHz and with an upgrade to 85 — 105 GHz 
operation in progress. Its location makes it a unique in- 
strument to study molecular-line sources in the southern 
hemisphere at high angular resolution, and we have re- 
cently used it to conduct a study of circumstellar HCN 
emission from the carbon star R Scl lii^ong ct al. 2o3)- 
The analysis in Sect. 3 and comparison with envelope 
models in Sect. 4 are carried out in the ww-plane in or- 
der to maximize the sensitivity and resolution of the data. 
The modelling is followed by a discussion in Sect. 5 and 
conclusions are presented in Sect. 6. 



2. Observations and data reduction 

At the time of the observations, the ATCA had three 
antennas of 22 m diameter equipped with dual polarisa- 
tion 3-mm receivers covering the bands 84.9 — 87.3 and 
88.5 - 91.3 GHz. We observed R Dor and Pup between 
August and October 2003 in three different array config- 
urations: EW367 (basehnes of 45, 90, and 135 m), 750B 
(60, 165, and 225 m), and H214 (75, 135, and 210 m). 
EW367 and 750B provided east-west baselines whereas 
H214 provided north-south baselines. The observations in 
each configuration occurred on single nights. The com- 
bined visibility plane coverage is shown in Fig. ^ 

All observations were conducted in clear to partly 
cloudy weather, with above-atmosphere single-sideband 
system temperatures of T*yg 350-500 K near the zenith. 
The correlator was configured to receive both linear polar- 
isations in two frequency windows: a narrowband (spec- 
tral line) window centred on the 'thermal' SiO v = 0, J = 
2^1 line at 86.847 GHz with 128 channels across 32 
MHz, and a wideband window centred on the 'maser' SiO 
?; = 1,J = 2^1 line at 86.243 GHz with 32 channels 
across 128 MHz. The pointing and phase centre was at 
aaooo = 04''36'"45!67, (Jaooo = -62°04'37'.'9 for R Dor and 
at aaooo = 07hl3'"32!31, ^2000 = -44°38'24'.'l for Pup. 
The field of view is determined by the ATCA primary 
beam which has a FWHM of about 36" at 86 GHz. 



^ The Australia Telescope is funded by the Commonwealth 
of Australia for operation as a National Facility managed by 
CSIRO. 
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Fig. 1. Coverage of the visibility plane for R Dor and 
Pup obtained after combining the three ATCA configura- 
tions used (EW367, 750B, and H214). 



The amplitude and phase gains were transferred from 
the maser to the spectral line window after applying scal- 
ing factors determined from a half-hour integration on the 
bandpass calibrator (3C279 or PKS B1921-293), which 
was also used to determine the channel-dependent gains. 
In the case of Pup, the transfer of the phase gains 
could not be performed successfully because of an inconsis- 
tency in the frequency setup between the source and band- 
pass calibrator; consequently we have applied phase self- 
calibration to the thermal SiO emission from this source 
assuming that it is circularly symmetric. For an interfer- 
ometer the astrometry is determined by the quality of the 
phase calibration. For R Dor the phase errors were ~ 5° 
and one would expect the position error to be about 0.1" 
for a 100 m baseline. 

Uranus was used to set the flux scale, assuming a uni- 
form disk with brightness temperature of 134 K. However, 
for the EW367 array, where the Uranus observation oc- 
curred in poor weather, the flux scale was tied to that of 
the 750B observation by assuming the SiO maser fluxes 
from both sources were unchanged over the 17-day gap 
between the observations. Assuming that any actual flux 
variations for the two sources were uncorrelated, the de- 
rived gains suggest an uncertainty in the flux scale for 
this configuration of ~ 30%. We also adjusted the antenna 
pointing once an hour on the SiO masers; typical pointing 
shifts were « 5". 

All d ata processing wa s conducted using the MIRIAD 
package l|Sault et alJl995|) . The calibrated visibilities were 
Fourier transformed using uniform weighting and a 0'.'5 
pixel size with a channel spacing of 1 kms""'^ (roughly 
the effective velocity resolution given the original channel 
spacing of .25 MHz pa 0.8 kms~^). The maps were then 
CLEANed l|Hogboml Il97^ down to a 2cr level over the 
inner 20"x20". 

The actual analysis and comparison with the model 
will be carried out in the uw-plane to maximize the sensi- 
tivity and resolution of the data. Thus we expect to obtain 
usable information on scales as low as 1", corresponding 
to the longest baselines. 



As an initial step in the reduction of the 750B and 
H214 data, an elevation-dependent gain curve was applied 
based on observations of an SiO maser taken on 2003 
Sep 4. No such gain curve was available for the earlier 
EW367 observation. Gain calibration as a function of time 
was then performed by assuming the SiO maser emission 
comes from a point source located at the phase centre. 
Since the maser emission is known t o emanate from a re - 
gion only a few AU in radius (e.g., ICotton et alJl200^ . 
this is a reasonable assumption for our resolution of a few 
arcseconds (1 AU — 20 milliarcseconds at the distance of 
our nearest source, R Dor). To derive the amplitude gains 
we also assumed that the integrated SiO maser flux was 
constant throughout each observation. We found that the 
variation in the gain was generally dominated by an un- 
corrected gain-elevation dependence, limiting any intrinsic 
flux variation to < 10%. 



3. Observational results 

3.1. R Dor 

The SiO w = 0,J=2^1 velocity channel maps of R 
Dor, binned into 2 kms^^ channels, are shown in Fig. |2| 
The synthesised beam is 2'.'6 x 2'.'0 with a position angle of 
—56° and the RMS noise in the 1 kms^^ velocity channel 
maps is (7 = 30 mJybeam"^, as estimated from regions 
free of source emission. The signal-to-noise ratio is high 
with the peak emission reaching a level of about 300(7. 

The velocity channel maps suggest that the SiO v = 
0, J = 2 — > 1 emission towards R Dor is moderately 
resolved, a result which is confirmed by the visibility 
analysis (see below). The brightness distribution appears 
to have an overall circular symmetry. From fitting two- 
dimensional Gaussian brightness distributions to the in- 
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RA offset (arcsec) RA offset (arcsec) RA offset (arcsec) 

Fig. 2. Velocity channel maps of SiO u=:0, J = 2^1 emission from R Dor. The contour levels are 0.03n Jy beam~^, 
where n = —9,-3,3,9,27,81,243 (negative values have dashed contours), and the beam size is 2'.'6 x 2'.'0 with a 
position angle of —56° as indicated in the lower right hand corner of each panel. The velocity channels (given in the 
LSR frame and indicated in the upper left corner) have been binned to 2 km s^^. The systemic velocity is 7 kms~^ 
as determined from CO observations. Offsets in position are relative to 0:2000 = 04'^36™45!67, S2000 = — 62°04'37'.'9. 



dividual channel maps, it is concluded that no systematic 
variations in position with velocity are present. The in- 
dividual offsets are smaller than O'.'l, within the pointing 
uncertainty (Sect.EJ, confirming the adopted source posi- 
tion. Across the peak of the line profile, from 6 — 8 km s^^ 
(LSR), the average values of the major and minor axes 
are 3'.'3 and 3'.'7 (FWHM), respectively. The beam has not 
been deconvolved from these values. 

As interferometers lacks sensitivity to large scale emis- 
sion it is of interest to find out whether part of the emission 
is resolved out in the ATCA data. In Fig. El the ATCA 
SiO data has been restored with a 57" circular beam 
(solid line; this represents the beam of the Swedish-ESO 
Submillimetre Telescope, SEST). A comparison with the 
equiva lent SEST spectrum from iGonzalez Deleado et "all 
l|2n0 3h (histogram) shows a striking resemblance both with 
regard to the peak flux (« 20 Jy) and the overall fine 



shape. This strongly suggests that all the flux is recov- 
ered by the ATCA interferometer, at least within the 57" 
SEST beam, and that virtually all of the SiO emission 
comes from the central component. The SEST spectrum 
has been converted from main beam brightness tempera- 
ture scale to Jy using 

S = ?7mbr~^T,„b, (1) 

where the main-beam efficiency rj^h = 0.75 and the sen- 
sitivity = 25 JyK~^. The telescope parameters are 
taken from the SEST homepage^. 

Assuming that the emission has an overall spherical 
symmetry, circularly symmetric Gaussians have been fit- 
ted fitted to the visibilities. Given the higher sensitiv- 
ity in the ww-plane there appears to be a systematic 

^ www . Is . eso . org/lasilla/Telescopes/SEST/SEST . html 



F. L. Schoier et al.: Probing the inner wind of AGB stars 



5 



20 



10 




10 



10 

[km s~^] 



20 



Fig. 3. SiO t; = 0, J = 2 "> 1 spectra for R Dor. 
The solid line is the SEST single dish observation 
ijGonzalez Deleado et al.l l2003(l . whereas the histogram 
shows the spectrum at the phase centre derived from the 
ATCA data by restoring with the SEST beam of 57". 
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Fig. 4. Offset from the phase centre of the velocity chan- 
nels as determined from fitting circular Gaussians to the 
R Dor data in the uw-plane. Indicated are the velocities 
of the two extreme channels (0 and 14 kms^^) for which 
reliable fits could be obtained. 



trend with velocity as shown in Fig. 0] The largest off- 
sets are measured at the extreme velocities, i.e., the blue- 
and redshifted emissio ns are displaced from each other. 
iLindavist et al.l l)200(j ). when analysing several molecular 
line emissions from the high mass loss rate carbon star CIT 
6, also noted that there was a common trend in position 
with velocity. A reasonable explanation for this behaviour 
is that the CSE is asymmetric to some degree, perhaps a 
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Fig. 5. Size (FWHM) of the observed SiO emission to- 
ward R Dor as function of velocity channel (given in 
LSR frame), averaged in 1 kms~^ bins, as estimated from 
Gaussian fits to the azimuthally averaged visibility data. 



bipolar outfiow, but the present angular resolution is not 
enough to make this apparent in the channel maps. Fig. 13 
The resulting FWHM of Gaussians fitted to the az- 
imuthally averaged visibilities are plotted in Fig. |S1 The 
overall variation in size with line-of-sight velocity is as ex- 
pected for a well resolved expanding envelope, where gas 
moving orthogonal to the line of sight subtends a larger 
solid angle than radially moving gas near the extreme 
velocities. The largest angular extent of the emission is 
3'.'0 ± 0'.'2 in the velocity range 6 - 8 km s^^ (LSR). This 
corresponds to a radial size of 1.0 x 10^^ cm at the dis- 
tance of 45 pc. We note that a single Gaussian provides 
a reasonable fit to the visibilities on all baselines except 
those close to the systemic velocity at 7 kms^^ (LSR). 
Introducing an offset in the flux scale, representing an un- 
resolved additional component, significantly improves the 
fit at those velocities. In Sect. 14.31 it is shown that the vis- 
iblilities obtained for this source are best fitted using two 
components: one compact with high a SiO abundance and 
one more extended with a significantly lower abundance. 

3.2. Pup 

The SiO w = 0, J = 2— +1 velocity channel maps of L^ 
Pup, binned into 2 kms^^ channels, are shown in Fig.|H| 
The synthesised beam is 2'.'5 x 1'.'9 with a position angle of 
14° and the RMS noise in the 1 kms~^ velocity channel 
maps is cr = 40 mJybeam"^, as estimated from regions 
free of source emission. The signal-to-noise ratio is high 
with the peak emission reaching a level of about 60cr. 

As expected from the single-di sh modelling performed 
by iGonzalez Delgado et alJ l)2003|) L^ Pup has a smaller 
envelope than R Dor. Since we had to self-calibrate the 
data in phase for this object, we have assumed that the 
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Fig. 6. Velocity channel maps of SiO u = 0, J = 2^1 emission from \? Pup. The contour levels are 0.04n Jy bearn"^, 
where n = —2,2,4,8,16,32 (negative values have dashed contours), and the beam size is 2'.'5 x 1'.'9 with a position 
angle of 14° as indicated in the lower right hand corner of each panel.. The velocity channels (given in the LSR frame 
and indicated in the upper left corner) have been binned to 2 km s~^. The systemic velocity is 33 kms~^ as determined 
from CO observations. Offsets in position are relative to q;2ooo = 07''13™32!31, (520oo = -44°38'24'.'l. 



source distribution is symmetric and centred on the stellar 
position. Across the peak of the line profile, from 33 to 35 
kms~^ (LSR), the average values of the major and minor 
axes are 2'.'9 and 2'.'2 [FWEM), respectively (note that the 
beam has not been deconvolved from these values). 

In Fig.Elthe ATCA SiO data for \? Pup has been re- 
stored with a 57" circular beam (solid line), representing 
the SEST be am. A comparison with the equ ivalent SEST 
spectrum of iGonzalez Deleado et al.l l)2003fl (histogram; 
converted from main beam brightness temperature scale 
to Jy using Eq. ^ shows a striking resemblance both with 
regard to the peak flux (« 3 Jy) and the overall line shape. 
As for R Dor this strongly suggests that all the flux is re- 
covered by the ATCA interferometer, and that all of the 
SiO emission comes from the central compact component. 

The result of fitting Gaussians to the azimuthally aver- 
aged visibilities is presented in Fig.|Hl The variation in size 



is also here consistent with a resolved expanding envelope. 
The largest angular extent of the emission is 1'.'6 ± O'.'l in 
the velocity range 33 — 35 kms~^ (LSR). This corresponds 
to a radial size of 1.0 x 10"'^^ cm at the distance of 85 pc. 
Contrary to R Dor, a Gaussian gives a reasonably good 
description of the visibilities on all baselines and for all 
velocity channels. 

4. Radiative transfer modelling 

A detailed non-LTE ra diative transf er code, based on the 
Monte Carlo method l|Berneslll979|) . is used to perform 
the excitation analysis and model the observed circum- 
stellar line emission. The code is described in detail in 
ISchoier fc Olofssonl l)200l|) and has been benchmarked, to 
high accuracy, against a wide varietv of molecular l ine ra- 
diative transfer codes in Ivan Zadelhoff eT 
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Fig. 7. SiO V = 0,J = 2 1 spectra for 

Pup. The solid line is the SES T single dish observation 
ijGonzalez Deleado et all l2003(l . whereas the histogram 
shows the spectrum at the phase centre derived from the 
ATCA data by restoring with the SEST beam of 57". 
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Fig. 8. Size {FWHA4) of the observed SiO emission to- 
ward Pup as function of velocity channel (given in 
LSR frame), averaged in 1 kms^^ bins, as estimated from 
Gaussian fits to the azimuthally averaged visibility data. 



4.1. The physical structure of the CSE 

The CSEs are assumed to be spherically symmetric and 
to expand at constant velocity, and they are produced by 
constant mass loss rates. Both R Dor and Pup were 
included in the large survey of semiregular M-type AGB- 
variables bv lKerschbaum fc Olofsson (1999.) . Their overall 
wind properties such as mass loss rate, terminal expan- 
sion velocity and temperature structure as constrained by 
single-dish CO millimetre line emission were presented in 



lOlofsson et alJ l)2002|) . where also a more detailed descrip- 
tio n of the circurn s tellar model is given. 

lOlofsson etliil l|2002l) noted that the CO lines for 
Pup are better reproduced when adopting a higher value 
of 1.0 kms^^ for the turbulent velocity than the nominal 
value of 0.5 kms~^. Further support for a larger turbu- 
lent velocity are the SiO line wings that are commonly ob- 
served in the low expansion velocity sources such as R Dor 
and Pup I Gonzalez Delgado et alJEflO,'^ . Therefore, in 
the present analysis the turbulent velocity width is used 
as a free paramete r . 

lOlofsson et all ll2002^ calculated the kinetic temper- 
ature distribution including, e.g., CO line cooling. Here 
we add cooling by H2O since this may be of impor- 
tance i n the region where the SiO l ine emission is pro- 
duced ( Gonzalez Delgado et al.lEoO.'^ . The importance of 



H2O as a coolant in t he inner wind was also noted by 
IZubko fc Elitzuil l)2000|) in their study of W Hya. In order 
to be consistent we have therefore produced new circum- 
stellar models for also the CO line emission. The inclusion 
of H2O line cooling lowers the temperature significantly 
(from 15 — 40%) in the inner part (r < 1 x 10^^ cm) of 
the envelope, but has only a little effect on the mass loss 
rate derived from CO observations. More details on how 
the H2O hue cooling is treated , and its effect on CO and 
SiO line intensities, are g i ven inlOlofsson et all l(2002l ) and 
iGonzalez Delgado et all l)2003(l . As an example, we note 
that excluding H2O line cooling increases the predicted 
SiO line intensities in our best fit model for R Dor (pre- 
sented in Sect. Ol by 10 - 15% for the J = 5 ^ 4 and 
J = 6 ^ 5 SEST observations as well as the J = 2 ^ 1 
interferometric ATCA observations. These are the data 
most sensitive to the physical and chemical conditions pre- 
vailing in the inner wind. For comparison, the predicted 
SEST J = 2 ^ 1 and J = 3 ^ 2 hnc intensities only 
increase by a few %. 

We find that the molecular line emission is better re- 
produced (in particular that of SiO) using a turbulent 
width of 1.1 kms~^ for Pup and 1.5 kms^^ for R Dor. 
Adopting a larger turbulent width means that the expan- 
sion velocity needs to be lowered somewhat in order to 
still provide good fits to the line widths. For for Pup 
an expansion velocity of 2.1 kms~^ was found and for 
R Dor 5.3 kms~^. These values are about 20% low er than 
those presented in ^Gonzalez Delgado et alJ l)2003l) where 
a much lower microturbulent velocity of 0.5 kms~^ was 
adopted. Changing the prescription of the velocity field, 
and in particular the turbulent velocity, changes the exci- 
tation conditions and somewhat different mass loss rates 
are obtained. For Pup a mass loss rate of 2.7 x 10^^ 
Mq yr~^ is found to provide a very good fit using an h- 
parameter of 0.05 typical of low mass loss rate objects. 
The /i-parameter contain much of the properties of the 
dust, such as size and density of individual grains and the 
dust-to-gas ratio, and enters in the heating term in the 
energy balance equation which is solved self-consistently 
in the excitation analysis (see Schoier & Olofsson 2001 
for a more elaborate explanation of the /i-parameter). For 
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Table 1. Source parameters obtained from modelling CO 
millimetre line emission. 
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From lQlofsson et al.l (|20o3) 
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Fig. 9. Azimuthally averaged visibility amplitudes, in 5 
kA bins, across the peak of the line profile (6 — 10 kms^^ 
in the LSR frame) for R Dor. The error bars indicate the 
2(T statistical error and does not contain the calibration 
uncertainty of about 15%. The solid line indicates the best 
fit model (Xq = 4 x 10"^ = 1.2 x lO^^ cm; = 
3 X 10^6; re = 3.3 x 10^^ cm), the da s hed li ne the best 
fit model from lGonzalez Delgado et aP l|2003|) {Xq = 5 x 
10~^; fe = 3.3 X 10"'^^ cm), and the dashed-dotted line the 
best fit model without freeze-out {Xq — 5.5 x 10~^; = 
1.0 X 10^5 cm). The dotted histo gram indicates the zero- 
expectation level, i.e., when no source emission is present. 



state. Relevant molecular d ata are summarized in 
iGonzalez Dekado et all l)2003l) . 

The abundance distribution of the circumstellar SiO is 
of importance in the radiative transfer modelling, and it is 
determined by two processes, dust condensation and pho- 
todissociation, as the gas expands away from the star. A s 
outlined in more detail in lGonzalez Delgado et ahl l)2003|) . 
the abundance is expected to decline, from the stellar at- 
mosphere equilibrium chemistry value, once condensation 
of SiO onto dust becomes efhcient when the kinetic tem- 
perature drops below the condensation temperature. This 
requires a minimum density which is not obtained in the 
low mass loss rate objects. For the higher mass loss rates 
the abundance declines drastically until the condensation 
ceases to be effective. Eventually, photodissociation sets 
in and this limits the size of the SiO envelope. Thus, to a 
first approximation the SiO abundance distribution con- 
sists of two components, a compact high abundance region 
close to the star, and a more extended low abundance re- 
gion. The former dominates for low mass loss rates and 
the lat ter for high mass loss rates. iGonzalez Delgado et alJ 
l|2003l) used only one component and its size was estimated 
from a formula obtained by fitting multi-line SiO data for 
eleven stars. 

In order to investigate which SiO abundance distribu- 
tion that best reproduces the flux picked up by the in- 
terferometer, the same (u, v) sampling was applied to the 
predicted brightness distribution from the model envelope. 
Two types of abundance profiles were used (by abundance 
we here mean the fractional abundance of SiO with respect 
to H2): i) a simple Gaussian abundance profile with an in- 
ner abundance and an e-folding distance and ii) a 
Gaussian profile with a depleted abundance Xy> and an 
e-folding distance Tc, and an added high abundance in- 
ner component Xq introduced as a step function inside a 
radius r^- This provides a first approximation to the sce- 
nario outlined above where a combination of adsorption 
of SiO onto dust grains and photodissociation dictate the 
abundance profile. 

4.3. R Dor 



R Dor a mass loss rate of 1.2 x 10 ^ Mq yr ^ is obtained 
with an /i-parameter of 0.05. 

The resulting envelope parameters are summarized in 
Table ^ The distances adopted are from lOlofsson et al.l 
l|2002(l and are based on the assumption that the stellar 
luminosity is 4000 Lq . For a discussion on this in relation 
to the measured Hipparcos distances, see lOlofsson et all 
ll2002l) . 

4.2. SiO model 

The excitation analysis includes radiative ex- 
citation through the first vibrationally excited 



The SiO abundan ce d istribution used by 
IGonzalez Delgado et all l|2003l) {Xq = 5 x 10"*^; 
Te = 3.3 X 10^^ cm) fails to explain the observed 
visibilities on all baselines, see the dashed line in Fig. 
This indicates that the abundance is higher in the inner 
part of the envelope than what is given by the model. 
It should however be pointed out that the analysis of 
IGonzalez Delgado et al is based on a size formula 

for the SiO envelope which is obtained through multi-line 
fitting of data for eleven stars. A better fit is obtained 
for Xo = 5.5 X 10~^ and re = 1.0 x 10^^ cm (dash-dotted 
line), but all the flux is not recovered on the shortest 
baselines using this model. Also, the flt to the single-dish 
data is worse. Instead, a combination of a high abun- 
dance compact component and a low abundance more 
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Fig. 10. Best fit model (solid line) overlayed on observations (histograms) in the form of single-dish SEST data 
ijGonzalez Deleado et al.ll2003tl as well as the spectrum at the centre pixel in the cleaned ATCA map for R Dor. The 
velocity scale is given in the LSR frame. The calibration uncertainty in the single-dish data is about ±20%. 
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Fig. 11. Azimuthally averaged visibility amplitudes, in 5 kA bins, across various parts of the line profile (as indicated 
in the upper right hand corner in kms~^ given in the LSR frame) for R Dor. The error bars indicate the 2(7 statistical 
error and does not contain the calibration uncertainty of about 15%. The solid line represents the best fit model using 
a constant expansion velocity. The dotted histogram indicates the zero-expectation level, i.e., when no source emission 
is present. 
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extended component {Xo = 4 x 10 



c = 1.2 X 10^5 cm; 

Xu = 3 X 10^^; To = 3.3 x 10^^ cm) provides a very good 
fit to the interferometer data (solid line), and at the same 
time to the single-dish data as seen in Fig. IIUI 

Of considerable interest is to investigate the assump- 
tion made in our models of a constant velocity throughout 
the SiO line emitting region. The observed line profiles 
are best fitted using a somewhat lower expansion veloc- 
ity of 4.9 kms~^ than obtained from the CO analysis (5.3 
kms""*^). This suggests that the wind has already been ac- 
celerated to within 5 — 10% of its terminal value at radial 
distances of « 1 x 10^^ cm 30 stellar radii). In Fig.lTTl 
the model predictions are compared with the observations 
in eight velocity intervals. The model can account for the 
trend that the gradient in flux with baseline decreases to- 
wards the line wings. To test the sensitivity of the SiO 
line emission to the adopted velocity law a velocity gradi- 
ent a ppropriate for a dust driven wind fe.g.. lHabing et all 



v{r) 



(2) 



was first assumed. The ratio of the terminal wind veloc- 
ity Voo = 5.3 kms~^ to the velocity vi at the inner ra- 
dius n = 1 X 10^^ cm was fixed to 0.25. This velocity 
law provides a good fit to the line width of all available 
SiO and CO single-dish spectra. It was found that this 
velocity law only changes the results very marginally and 
provides an equally good fit to the observed visibilities as 
does a constant velocity field. Models fo r pulsation driven 
winds l|Winters et alJ2000blEooll2003^ generally suggest 
a more slowly increasing velocity with radius. A velocity 
law having a constant velocity of 5.3 kms~^ outside 30 
stellar radii (« 1.4 x 10^^ cm) and 



v{r) = 5.3 



1.4 X IQiScm 



0.5 



kms 



(3) 



inside this region approximates such a slow increase. It 
provides equally good fits to the observations. The reason 
that no significant difference is observed is related to the 
fact that the SiO ti = 0, J = 2^1 emission is optically 
thick. Thus, the emission emanates from a small region 
close to Tc over which no large velocity gradient exist. 

A point worth mentioning is that the relatively high 
turbulent velocity of 1.5 kms^^ obtained for R Dor, as 
required to fit the line shapes for both the CO and SiO 
line emission, also significantly affects the excitation of, in 
particular, the lower rotational levels. In the modelling it 
was found that the norn inal value of 0.5 kms^^ used by 
IConzalez Delgado et al.l (|2003) results in population in- 
versions in the lower level transitions of SiO, i.e., maser 
action. The line profiles in that case have strong peaks at 
the extreme velocities which is not observed. The higher 
turbulent velocity helps to quench the maser emission and 
produces line profiles that very closely resembles the ob- 
served ones, Fig. EH 
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Fig. 12. Averaged visibility amplitudes, in 5 kA bins, 
across the peak of the line profile (32 — 36 kms^^ in the 
LSR frame) for Pup. The error bars indicate the 2a 
statistical error and does not contain the calibration un- 
certainty of about 15%. The sohd fine indicates the best fit 
model (Xo = 6 X 10"^ rc = 8 x lO" cm; Xd = 2 x 10"^ 
Te = 2.1 X 10^^ cm), the das h ed lin e the best fit model 
from iGonzalez Delgado et all l|2003|) (Xq = 1.4 x lO'^; 
To = 2.1 X 10^^ cm), and the dashed-dotted line the best fit 
model without freeze-out (Xq = 5 x 10~^; — 1.1 x 10^^ 
cm). The dotted histogram indicates the zero-expectation 
level, i.e., when no source emission is present. 



4.4. /.2 Pup 

Pup has a significantly lower (a factor of « 4) mass 
loss rate than R Dor. However, the expansion velocity of 
the wind is also a factor of two lower resulting in only a 
factor of about two lower density in the wind. Thus, the 
SiO envelope characteristics should be rather similar for 
the two sources. 

The SiO abund a nce distribution used by 
iGonzalez Delgado et all ll2003l) (Xq = 1.4 x 10"^ 
Te = 2.1 X 10^^ cm) does not provide a good fit to the 
observed visibilities as illustrated by the dashed line 
in Fig. El While the model can account for the total 
fiux, as indicated by the shortest baselines, it fails to 
account for all the flux at the longer baselines. This again 
indicates that the abundance is higher in the inner part 
of the envelope than what is given by the model. The 
flux seen by the interferometer at the longer baselines 
can be accounted for by raising Xq to 4.5 x 10^^, at 
the same time that the size of the region needs to be 
lowered to 1.1 x 10^^ cm in order to reproduce the total 
flux seen at the shorter baselines. This model is shown in 
Fig. El as the dashed-dotted line and account very well 
for the flux at all baselines. It should be noted that this 
model provides a somewhat worse fit to the single-dish 
data, with a reduced Xrcd '^^ ^-^^ particular the high- J 
lines come out too strong in the model by about 50% on 
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Fig. 13. Best fit model (soli d line) overlayed on observations (histograms) in the form of single-dish SEST 
ijGonzalez Delgado et al.l|2003|) and JCMT data (unpublished data) as well as the spectrum at the centre pixel in 
the cleaned ATCA map for Pup. The velocity scale is given in the LSR frame. The calibration uncertainty in the 
single-dish data is about ±20%. 
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Fig. 14. Averaged visibility amplitudes, in 5 kA bins, across various parts of the line profile (as indicated in the upper 
right hand corner in kms~^ given in the LSR frame) for Pup. The error bars indicate the 2a statistical error and 
does not contain the calibration uncertainty of about 15%. The solid line represents the best fit model using a constant 
expansion velocity. The dotted histogram indicates the zero-expectation level, i.e., when no source emission is present. 



average. In the light of the results obtained for R Dor, 
the combination of a high abundance compact compo- 
nent and a lower abundance more extended component 
(Xq = 6 X 10-^ Tc = 8 X 10" cm; Xd = 2 x 10"^; 
fe = 2.1 X 10^^ cm) produces an equally good fit to the 
interferometer data (Fig. solid line) and at the same 
time provides a much better fit to the available single-dish 
data (x^ed = 1-5)' Fig. El We note that our model can 
not reproduce the observed J = 5 ^ 4 line at the SEST 
which comes out with an integrated intensity of about 



30% higher than observed. The absolute calibration at 
the SEST is expected to be better than approximately 
±20%. 

In Fig. 1141 the model predictions are compared with 
the observations in four velocity intervals. The model can 
account well for the trend that the gradient in flux with 
baseline decreases towards the line wings. Note that the 
absolute scale of the flux in the model is somewhat lower in 
the wings. This is due to the fact that the model does not 
reproduce the line wings perfectly. Although, the source 
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size estimated from the line wings is very well represented 
by the model. As for R Dor there is no evidence of any 
large gradients in the velocity field. Note that also for 
Pup the SiO i) = 0, J = 2^1is optically thick and hence 
can not be expected to contain detailed information of the 
velocity law on smaller scales. 

5. Discussion 

5.1. SiO abundance distribution 

Our modelling of the interferometer SiO data clearly show 
the need of an abundance distribution with a high abun- 
dance {Xq w 4 X 10~^), compact (« 1 x 10^''^ cm) compo- 
nent in order to reproduce the visibilities obtained with 
ATCA for both Pup and R Dor. In addition, there is 
evidence of a low abundance (Xd ~ 2 — 3 x 10^^), more 
extended (w 2 — 3 x 10^^ cm) component. 

This scenario is very similar to the interpretation 
of interferometric S iS observations of the carbon star 
IRC-l-10216 made bv lBieging k Tafallal l)l993(l . where the 
abundance of SiS is seen to drastically drop by two orders 
of magnitudes at a radial distance of about 2 x 10^^ cm. 

The values of the initial SiO abundance derived 
here is very close to that expected from stellar atmo- 
sphere chemistry in thermal equilibrium, 3 — 4 x 10~^ 
l|Willacv fc MillaJll997HDuari et al.lll999|) . Assuming so- 
lar abundances a maxi mum SiO fractional abundanc e 
of 7 X 10"^ is obt ained ijGonzalez Dekado et alJl200^ . 
iDuari et al ] lll999 h further showed that non-equilibrium 
chemistry, where the chemistry is controlled by shocks, 
only has a very limited effect on the abundance of SiO 
in the inner part of the wind. Thus, it is expected from 
theory that the photospheric abundance of SiO remain 
relatively unaffected until the combined effect of adsorp- 
tion and photodissociation sets in. Our data support such 
a scenario. This will be discussed further in Sect. 15.21 

As an alternative explanation to the low abun- 
dances of c ircumste llar SiO t hat we re estimat ed by 
iLambert fc~vanden Bou^ lll978|) and iMorris fc Alcockl 
l)l977|) . IScalo fc Slavskvl l)l98(]() suggested that photodis- 
sociation by chromospheric radiation may be important. 
This is clearly not the case for the two sources studied 
here but a statistical sample needs to be investigated in 
order to rule out such a scenario. Unfortunately, not much 
is known about the chromospheres of AGB stars. 

5.2. Grain formation 

The presence of a high initial SiO abundance followed 
by a drastic decrease indicates that adsorption onto dust 
grains could pl ay an important role in removin g SiO from 
the gas phase. iGonzalez Delgado et al.l l)2003l) presented 
a relatively simple model to estimate how effectively SiO 
molecules freeze-out onto dust grains. Using this we find 
condensation radii of « 5 x 10^^ cm for both stars (es- 
sentially determined by where the dynamical time scale 
equals the evaporation time). The degree of adsorption 



critically depends on the amount of dust available. For the 
dust parameters hinted at from the CO modelling and con- 
tained in the ^.-parameter we expect the SiO abundance 
to drop by about a factor of ~ 3 beyond this radius (for 
both stars). Likewise, we can use the simple photodisso- 
ciation model outlined in lGonzalez Delgado et aP l|2003fl 
to derive SiO photodissociation radii of about 10^^ and 
3 x 10^^ cm for Pup and R Dor, respectively. Thus, 
within the considerable uncertainties, we expect the SiO 
line emission to be strongly dominated by the compact 
high abundance region in Pup. In the case of R Dor the 
low abundance region between the condensation and pho- 
todissociation radii will also contribute. Our observational 
results are certainly consistent with this, even though the 
estimated abundance drop beyond the condensation ra- 
dius is larger than predicted by the model. Given the sim- 
plicity of the treatment of the adsorption process and the 
uncertainty in the adopted dust parameters, we find this 
acceptable. More elaborate chemical models than the ones 
presented bv lGonzalez Delgado et al.l ( 2003? ) are needed to 
take the combined effect of adsorption and photodissoci- 
ation properly into account. 

In addition, as indicated in Sect. 15. li the role of chromo- 
spheric radiation as a means of dissociating S i O in these 
stars needs to be investigated. IScalo fc Slavskvl 1 19H(t sug- 
gested that SiS might be the best probe of the degree of 
adsorption of Si bearing molecules since it can be read- 
ily formed though reactions with H and H2 in the in- 
ner wind. In contrast there is no easy means of forming 
SiO. However, this was based on the assumption that all 
Si-molecules were destroyed by chromospheric radiation. 
Nevertheless, observations of SiS is an important comple- 
ment to the SiO ob servations, and it has been detected 
in 0-rich envelopes l)Lindavist et al.lll988H01ofsson et all 
li 



5.3. Wind dynamics 

The success of the modelling using a constant velocity 
(and a higher micro-turbulent velocity width) strongly 
suggests that the wind has reached its terminal velocity 
already within ~ 20 — 30 stellar radii. This is further sup- 
ported by the fact that the SiO and CO line data can 
be modelled using the same expansion velocity (within 
5 — 10%). The marked difference in shapes between the 
SiO and the CO lines is attributed to strong self absorp- 
tion of emission in the part of the wind moving towards 
the observer in the case of SiO. Th us, we support the 
conclusion bylSahai fc Biegind l)l993(l rather than that of 
iLucas et all l|l992j) where an extended region of acceler- 
ation was advocated, but also note that we have so far 
studied only two low mass loss rate objects. 

It has been su ggested from hydrodynamical calcula- 
tions l)Winters et al. 2000b. 2002, 2003) that in very low 
mass loss rate (M < 5 x 10"'' M0yr~^), low expansion 
velocity {< 5 kms~^) objects, such as R Dor and Pup, 
the main driving mechanism behind these tenuous winds 
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Table 2. Predicted SiO line intensities (integrated over 
the line and in main-beam brightness scale using a 12 m 
telescope) towards R Dor for selected rotational transi- 
tions. 
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65 


63 


30 -+ 29 


1300.5 
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52 
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Best fit model from lOonzalez Deleado et all (|2003h (Xq = 
5 X 10~^ rc = 3.3 X 10" cm). 

^ High abundance compact component (Xo — 5.5 x 10~^; ~ 
1.0 X 10" cm). 

Best fit model (Xo = 4 x 10"^ rc = 1.2 x lO" cm; Xd = 
3 X 10"^ rc = 3.3 X 10" cm). 



is stellar pulsation, and the wind acceleration is low. This 
is in contrast to winds which reach higher mass loss rates 
and expansion velocities where dust plays a primary role, 
and where the wind acceleration is very efficient. From the 
present data we can not distinguish between these com- 
peting scenarios mainly due to the SiO emission being op- 
tically thick, and therefore most of the emission emanates 
from a region where the wind is close to the terminal ve- 
locity. Observations of optically thin SiS lines could pos- 
sibly be useful to probe the velocity field. However, the 
SiS line emission is expected to be more than an order of 
magnitude weaker than that of SiO. 

5.4. Predictions for high- J SiO lines 

Presented here in Table El are estimated line intensities 
for selected transitions of SiO based on our models for 
R Dor using a 12 m telescope. Th e differences between the 
one-co mponent model derived bv lGonzalez Deleado et all 
l)200,ll) . and the new models required to fit the interferome- 
ter data presented here, increase as higher frequency tran- 
sitions are considered. This shows the potential of single- 
dish telescopes such as CSO, JCMT and the upcoming 
APEX"^ and ASTE telescopes to study the properties of 
stellar winds close to the central star. 

APEX and ASTE are of particular interest since they 
will eventually open up the THz region for observations. 

^ The Atacama Pathfinder Experiment (APEX) , is a collab- 
oration between Max Planck Institut fur Radioastronomie (in 
collaboration with Astronomisches Institut Ruhr Universitat 
Bochum), Onsala Space Observatory and the European 
Southern Observatory (ESO) to construct and operate a mod- 
ified ALMA prototype antenna as a single dish on the high 
altitude site of Llano Chajnantor. 



They will also be the only telescopes at the southern hemi- 
sphere being able to observe at sub-milimetre wavelengths. 

6. Conclusions 

High resolution, interferometric, millimetre observations 
at w 1" resolution of SiO v = 0, J = 2-^l line emission 
towards the two 0-rich AGB stars Pup and R Dor have 
been performed. The emission is resolved, very centrally 
peaked, and suggests an overall spherical symmetry, even 
though there is an indication of a departure from spheri- 
cal symmetry in the case of R Dor at the arcsecond scale. 
A detailed excitation analysis was performed suggesting 
that the SiO abundance is very high (« 4 x 10^^) in the 
inner part (^ 1 x 10^^ cm) of the circumstellar envelopes 
around both stars, consistent with predictions from LTE 
stellar atmosphere chemistry. For R Dor the interferome- 
ter data further suggests that there is an additional, more 
extended, component with a significantly lower SiO abun- 
dance. We interpret this as the result of effective adsorp- 
tion of SiO onto dust grains. Such a conclusion is less clear 
in the case of Pup. 

A comparison of model and observed line profiles fur- 
ther suggests that micro-turbulent motions are of the or- 
der of 1 — 1.5 kms^^. This is a significant fraction of 
the wind velocity in these slowly expanding winds (< 5 
kms^^), and it produces strong self absorption of the blue- 
shifted emission giving the SiO line profiles their charac- 
teristic shape. Additionally, the interferometer data pro- 
vide constraints on the size of the region in which the wind 
is being accelerated. It is found that a constant velocity 
over the SiO emitting region can successfully account for 
the visibility amplitudes obtained at different velocity in- 
tervals. Furthermore, it is possible to model the line shapes 
of both SiO and CO emission using the same expansion 
velocity (within 5 — 10%) for both our sources. This con- 
strains the acceleration region to within < 20 — 30 stellar 
radii. 

We conclude that SiO line emission plays a very im- 
portant role in the study of circumstellar envelopes, both 
with respect to its dynamics and the gas and grain chem- 
istry. However, it should be remembered that both sources 
studied here have very low mass loss rates and expansion 
velocities and that a more representative sample of 0-rich 
sources needs to be studied in order to solidify the conclu- 
sion reached here. This could also be used to investigate 
the role of chromospheric radiation in photodissociation of 
SiO molecules. Such an effort is currently underway, and 
it justifies the development of a more elaborate model to 
describe the abundance distribution of SiO. 
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